Hedging Swing contract on gas market 
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1 Introduction 

Swing options on the gas market are american style option where daily quantities 
exercices are constrained and global quantities exerciced each year constrained too. 
The option holder has to decide each day how much he consumes of the quantities 
satisfying the constraints and tries to use a strategy in order to maximize its expected 
profit. Based on the indexation principle, the pay off fonction is a spread between 
the spot gas market and the value of an index composed of the past average of some 
commodities' spot or future prices. Most of the time the commodities involved are 
gas and oil. 

When the index is a fixed one (then it is a strike), the valuation of Swing option 
is a classical problem solved by dynamic programming. In order to implement the 
dynamic programming method , first trinomial trees were used lfl3ll to estimate con- 
ditional expectation, then Longstaff Schwartz method [ 15 1 using regression and Par- 
tial Differential Equation methods [ 18 1. 

At the opposite, Swing options on the gas market are very high dimensional prob- 
lems due to the index definition. This problem can be linked to the problem of pric- 
ing moving average american options. These moving averaging options with early 
exercice features have been studing in [4|, [8|, [10]. A common approach (see |4|) 
is to use least square Monte Carlo to estimate conditional expectation in Monte- 
carlo algorithm. Recently [2| derived a new methodology based on exponentially 
weighted Laguerre polynomials expansion to approximate the moving average pro- 
cesses. The methodology developped is efficient but numerical results showed that 
the usual regression method regressing on both the gas price and the index value is 
accurate enough. 

If the question of the valuation of these contracts is rather well explored, the hedge 
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of such a contract has never been studied. Most of the time practionners just hedge 
the gaz spot component in the index ignoring the stochasticity of the index. Be- 
sides when practionners want to simulate the portfolio risk they have to include the 
dynamic hedge they plan for next years in order to assess their risk accurately. As 
for gas storage, a methodology based on tangent process has been developped by 
ifTTl to evaluate the Delta and follow the hedging strategy. The mathematical proof 
of this approach based on the Danskin's theorem is given in 0. In this paper the 
methodology based on the tangent forward method calculated during a dynamic pro- 
gramming algorithm is developped for the gas Swing options. Its accuracy is tested 
on a realistic contract where: 

• prices follow a two factor model (see for example lfT6ll ). 

• hedging products available are day ahead products till the end of the week, week 
ahead products till the end of the month, and monthly contract as they can be 
found at Henry Hub. 

The efficiency of the hedge and the accuracy of the calculations are studied depend- 
ing on : 

• the regressor used to calculate conditional expectation, 

• the products involved in the index that we decide to hedge, 

• the frequency of the hedge with respect to the products involved in the index 
including exchange rates. 

In the sequel we suppose that conditional expectation are calculated by the Longstaff- 

Schwartz method lfT31l adapted with local functions as explained in (6). 

In the next section[2] the contract and the price model are described. In the following 

section[3] the methodology used to simulate the dynamic hedging is explained. Next 

numerical results are given and a concluding section|5]gives some recommandations 

in particular in terms of frequency of the dynamic hedging and the components to 

hedge. 



2 The contract and the price modelization 

2.1 The contract 

The Swing contract is first characterized by the quantities that can be purshased. 
Daily exercice dates are given in days from t start to t enc i. Each day t the option holder 
can exercise a quantity q, satisfying 

< q t < q,„ax 

Moreover the global quantity consummed Q t is constrained 

Qmin — Qt — Qtnax- 

tstart<t<t en{ i 
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The unitary pay off (for a quantity one) is given at a day t by 

where S® is the gas price and I, an index constant on each interval [7], Tj + \ [ where 7} 
(z = to n) are update dates defined in the contract. These dates correspond to the 
beginning of some months and typical values for 7) are first day of each month, or 
first day of every two or three months. 

For each exercice day i, we note k(i) = max{y,7) < ;}. In the case of an additive 
index, the structure of the index for each exercice day i is defined by the affine sum 
of N components by : 

Ii = h m = ao + L a j (SJ T - bj)X]. 

;'=i 

where : 

• Sjj- — ij2* =1 5y j is the average of commodity j (potentially a forward 
price) on the m months expressed in days preceding the lag corresponding to 
^ months expressed in days involved in the index valid from date T (see figure[T|), 

• Xj is the exchange rate between the domestic currency and the foreign currency 
of the commodity j at date 7^, 

• aj, j = 0...N, bj = 1 , ...7V are coefficients given by the contract. 

Remark 1. In some simple cases the lag and averaging windows are the same for 
each component of the index and the contract windows can be decribed by the triplet 
m,l,p where p gives in months the validity period for the index, I the lag periode 
and m the averaging window length. 

Remark 2. In the formula, the exchange rates involved could be imposed by the 
contract at date T(k) — l(k) or an average of the exchange rate could be imposed. 



lag l(k) in months 



Averaging period m(k) in months 

*l(k) 0t(k) 

Fig. 1 Lag and averaging windows for calculation 
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2.2 Price Models 



All over this section, we shall consider a n-dimensional Brownian motion jz 1 , . . . ,z"} 
with correlation matrix p on a probability space (Q , J^, P) endowed with the natural 
(completed and right-continuous) filtration F = {^ t )t<T generated by {z 1 ,... ,z"} 
up to some fixed time horizon T > 0. The n value will depend on the different mod- 
els used for commodity prices and exchange rates. 



2.2.1 Future Price Model 

We suppose that the uncertainties in the commodity prices Sj, j — Q,N follows under 
the risk neutral measure a n-dimensional Omstein-Uhlenbeck process. The follow- 
ing SDE describes our uncertainty model for the forward curve F J (t, T) giving the 
prices of a unitary amount of a commodity j at day t for delivery at date T ( see Q): 

dF !p T ) = y a j{t)e- a ^-t)dz? +i , (1) 
FJ(t,T) ftx ' 

with a{ some volatility parameters and aj mean reverting parameters for commodity 
j- 

Remark 3. Most of the time a two factors model is used. In this model, the first 
Brownian motion describes swift changes in the future curve, the second one de- 
scribes structural changes in the gas market and deals with long term changes in the 
curve. The mean reverting parameter is generally taken equal to zero for the second 
term. 

With the following notations: 

z= l/l=i+l J 

W}' j = fakuSe-l^dg*', fori=l,...,n 
Jo 

the integration of equation ([T]) gives : 



-^Vj(t,T) + £e-^W;> J 



F J (t,T)=F J (t ,T)e 2 i=l . (2) 



With this modelization, the spot price for commodity j is defined as the limit of the 
future price : 
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S{ =]imF j (t,T). (3) 

As for the instantaneous exchange rate for commodity j, X, J , we choose a simple 
model 



jyj 

^ = MO - A(t))dt + a^dzf +J , (4) 
for j = l,..,N where 

• rl(t) is the short rate for the currency of commodity j, 

• r d {t) is the domestic short rate, 

• a J x is the volatility of the foreign currency associated to commodity j in domestic 
numeraire. 

In the sequel we suppose that the domestic interest rate ra(t) is for simplicity. 



2.3 Objective function 

Under the martingale probability, the owner of the option will try to optimize it's 
profits by maximizing the following 



J(q) = E 



tend 

E {*(s?-/,)} 

t—tstart 



where q = {q t ,t € [t s tart, ■■■,t en d]}> under the constraints : 

( q t G [0,qmax] 

l_ Qmin ^ Ldt=t start 1' — Qmax 

We introduce 

Qt= I qt 

t=tstart 

^ = {(%start> ■■■i ( lt end ),% G =^f,0 <q t < q m ax,Qmin < 2w — Qmax}, 

Z, = {zf,t < tj < n = 2(N+ l)n+Ar} 



and 



(5) 



(6) 



the stochastic part of the state vecteur at date f . 
The solution of the problem is given by J* : 

/*=max7( 9 ). (7) 

According to the Bellman principle, the optimal function value at date t expressed 
as a function of Q t and Z t satisfy : 
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r(t,Z,Q t )= sup {q{S°-I t )+E[J*(t + l,Z, + u Q r +q) \ Z, 

0<q< q max , 
Qt + ^Ltt>t Qmax ^ Qmin > 
Qt + q < Qmax 

(8) 

This approach is a classical one which has been used for example recently in [ 18 1. 



3 Optimization and dynamic hedging 

In the first part we give the equation for the dynamic hedge of each commodity and 
echange rate of the problem according to the methodology developped in [17|. We 
then explain why Delta have to be approximated and aggregated to be calculted on 
realistic problems. 



3.1 Continuous version of the hedge 

During the dynamic programming resolution, the daily quantities to exercice can 
be discretized if necessary or a bang bang approximation (which can be exact 
can be used. We note q*(t,Z t ,Qt) the optimal volume exerciced at date t when the 
quantity already consumed is Q, . We introduce Q*-' (z, c) the optimal volume level at 
date ii starting at level c at date i following a trajectory Z, with Z,- = z- The optimal 
volume is J^-mesurable and follows 

Q*i'(z,c) = C,yii < t s ta r t, ii > i, 

»-l 

Q*i\z,c) = c+ £ q*(k 1 Z k ,QY{z,c)) for t sturt < ii < t md . (9) 

k=tstartVi 

Thus Ylk=t , , 1*(k:Zk, Ql'°{z,0 j) corresponds to the sum of the optimal volumes 
exercised following the optimal strategy starting with a zero consummed volume at 
date with Zq = z- 

Following the methodology in lTT7ll we introduce the forward tangent process for 
commodity /' noted Y/' T satisfying : 

Y W=e-* vJ M + n,i>^'Wj = o t .„ )Ni (10) 
and for the exchange rate j the classical tangent process [9] : 

As prooved in [5 1, introducing the J^„, random variable : 
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D°(i,m,z,Q) = q*(m,Z m ,Q*J(z,Q))Y% m . (12) 

the sensibility of the contract value with respect to the arbitrage market (spot gas) 
for delivery at date m is then given by 

A°(i,m,z,Q) = E[D°(i,m,z,Q) | Z = z]/Y?< n ,t md > m > i (13) 
For the component j of the index, introducing 

&(i,m,l,z,Q) =q*(l,Z l ,Q;' i (z,Q))YA m , (14) 

and for m > i 

DJ(i, m ,z,Q)=ajZ £ P/(f,m,/,t,fl)^ w lrW "^^ ) ^ » 
* le[T(k)j{k+i)[ m W 

the sensibility of component /' of the index at date i and delivery at date m (m > i) 
is then : 

Aj{i,m,z,Q) = -E[D-'(i,m,z, Q) | Z, = z]/Y/' m (16) 



To ease the comprehension of formula ( 16 1, we explain it for the first index compo 



t=i t= m 

TTO-l(0)~m(fl3 i T(0)-\(0) T(2)-l(5-mf2J TPHK! 

T (°> pg^gB^psg^s^ | 

T(lJ-l(l>rn(lJ T(1>-1(U 

Fig. 2 Example of sensibility for first component of the index 

nent (commodity one) on figure [2] when the index is only updated tree times at date 
7b, T\ and T%, At date i, the first component of the index corresponding to commod- 
ity 1 with delivery at date m has a direct impact on the indexed price at which the 
gas will be paid for all the dates between 7b and T2 so the sensibility can be writ- 
ten as the sum of a component for delivery during the first and the second period. 

£ <7* Z/ , Q* [ ' J (z, Qj) for example is the optimal consumption on first peri- 
(e[r(o),r(i)[ 

ode exerciced with the index 7j-( ) constant on first period and thus the sensibility in 



equation ( 16 1 becomes : 



A\i,m,z,Q) = -^ UUE[ £ q\l,Z l ,Q]\z : Qmm Y ^ m \ Z i=^ 



le[T{0),T{\)[ 
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1 E[ £ q*(l,Zi,Q*'\z, Q))x£ m Y^ m | Z ; = z ] 
/e[r(i),r(2)[ / 

As for the exchange rate corresponding to component j of the index, introducing 

D J x (i,m,l,z,Q)=q*(l,Z l ,Q* l ' i (z,Q))Yl 

and 

D J x (i,z,Q)= aj £ L D x (iJ(k)J,z,Q)S^ (k \ 

k,T(k)>i le[T(k),T{k+l)[ 

the sensibility is given by 

A±(i,z,Q) = -E[D J x (i,z,Q) | Z, = Z ]/Y/ (17) 

Here for each period k that has not began, we assess the optimal volume exerciced on 
the period multiplied by the tangent process at T(k) for the exchange rate multiplied 
by the aj coefficient and the average value of the commodity j used for the index 
on period [T(k), T(k+ 1)[. The results are obtained by summing on the periods and 
taking conditional expectation. 



3.2 Discretized algorithm 

In order to solve equation rfSj, the methodology in Jfjl derived from the Longstaff- 
Schwartz methodology [ 15 1 can be used as in the case of American options. Using 
this methodology adapted to Swing option, optimal exercices and corresponding 
hedging strategies calculated at a date t depend on quantities consumed before this 
date. So optimal values and hedge have to be calculated on a grid used to discretize 
the possible volume consumed. Interpolation between Bellman values and hedging 
strategies is then used during the dynamic programming process. 
The volume is discretized on the grid 

Qi = lA, l = 0,...,L = Q max /A 

where A is the mesh size. Similarly to the case of American option, a Monte Carlo 
method is used to get some prices simulations (S/)-' for j = l,...,M at day ; = 
to t em i, for j = 0, ...,N. In the sequel (.) ; will stand for the simulation j of random 
variable (.). The conditional expectation E[ | Z,] has to be calculated in very high 
dimension but as shown in (2), it is possible to approximate the operator E[ | Z ( ] by 
the operator E[ | Sf ,/,•] with accuracy leading to regression in dimension 2. In the 
numerical part of the article a small study on the regressor will be achieved more 
thoroughly. In the sequel we note E[ | Z t ] an approximation of the conditional ex- 
pectation operator E[ | Z,] obtained by regression with |6| methodology where the 
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information contained in the stochastic state vector Z, has been approximated by 
some information contained in Z, vector. 

We note q* and Q* the estimation of the optimal volume q* and optimal volume 

levels Q* obtained by the Longstaff-Schwartz method. We note D°, L)i,b , D J X , 

D x , j = l,N the estimations of the variable D°, D^,D^, D J X and D J X on Montecarlo 
simulation for optimal control calculated q* . 

As explained in ifTTl . knowing D°(i+ l,m, (Z + i) J ,2/) f° r m = i+ l,...t eru j — 1, 
/ = 1 , . . . L, j = 1 , . . . M and using the equality 

Gti(&) J ,Q) - e*' i+1 ((^ + i) J ,G+ro',(^) J ,G)) 

for tend >m> i, the following backward recursion can be used to calculate D°(i, j, Zf,Q) 
at day i for all j > i, for all Q € {Q,/l =0,...,L}: 

{ D(i, i, (Z,) J \ Q) = q* (i, (Z)', Q) {Y?' j y 

\ £>{i,m,(Zi)},Q) = D(i+ l,m,(Z i+l )i ,Q + q*(i,(Zi)i ,Q)),m = i+lj end 

As for as the current backward recursion on option value, an interpolation on D 
values is necessary to estimate D(i + 1 , ., .) for the stock point Q + q* (i, (Z,) ; , Q). 
The Longstaff-Schwartz estimator of the conditional Delta is then evaluated for m > 
kby 

£°(i,m,ft)J,Q) = E[D(i,m,Zi,Q) | Z = (Z^}/ (Y°-"'y . (19) 

Similar recursion can be used for D J with j > 0. We first give the D J approxima- 
tion of D J : 

£>j(i,m,l,(Ziy\Q) =q*(l, {2tf,Q]\(ti)lQ))(Ylrf,i <m<l< t end (20) 
Introducing D J the D ; approximation : 



* ie[T(k),T(k+i)[ 



m(k) 



we then calculate the D values during the recursion easily supposing that D 3 (i- 
1,7, (Z;) J ,Q) have been previously calculated with 



I, 



D J (i, m, (Zi)J, Q) =D J (i+l, m, (Z i+I y, Q + q* (i, {%)}) ,Q),i<m 

(21) 

Equation (21 1 allows to calculate the hedge : 
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Aj(i,m,{Ziy,Q) = -E[D J (i,m,l,2 h Q) I % = (Z,) J }/(Y, j ' m ) J (22) 
As for the exchange rate we introduce D J X the D J X approximation : 

£> x {i,m,l,{2iy,Q) = q*(l, (Z,y,Q*- i ((Z i y\Q))(Y,^\i <m<l (23) 
and the D x approximation of D x 

b x (i,(Ziy,Q)=aj £ £ /^(i,r(*),/,(^) / ,G)(^)* ( V 

k/T(k)>i l€[T(k),T{k+l)[ 

(24) 

The following recursion adding at each day all the contribution on the sensibility 
due to the following day can be used during the backward recursion : 

' D x (i,(2 i )?,Q)=0,j>T(n), 

< b x (i,(Zi)J,Q) = ajZkSr(k)(i+l)Lie[T(k),nk+i)[Di(i,T(^^ 

+D X (i+l, (Z i+I y, Q + q* (i, (Z,y, Q)),i<T(n) 

(25) 

leading to the Delta approximation 

A x (i, fry, Q) = -E[D J x (i,Z„ Q) | Z, = (Z,-) W/) 7 (26) 

Using equatio n (fT~8] >, ( p~9] > sensibility with respect to gas is thus calcultated, using 
equation pT| ) \22\ sensibility with respect to commodity index components calcu- 
lated and using equation ( 25 1 and ( 26 1 sensibility to exchange rate calculated. 



3.3 Aggregation 

Not all product are available on the market. We note £i l , the set of future products 
available at date i for commodity /', and for all p € =2/, &p the delivery period 
associated to product p available for commodity j, r\ J p the beginning of the delivery 
period. Supposing that Vf > 0,V/? € =2/,Vj > i there exist £> j > p C J2l such that 

&p = Up e _ S j,p SP'p then it is possible to aggregate an approximated conditional Delta 
at date i per product with an ad hoc rule so that a dynamic programming approach 
is still usable. A first way to get the Delta on a delivery period p is to average the 



Deltas calculated with equations (19i and (22i calculated during recursion on the 
delivery period. This approach is memory cumbersome as noted in ifTTl . It is far 
more efficient to calculated the Delta directly on the product p during the backward 
recursion. [ 17 1 proposed to use in the continuous framework : 
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D°{i,p,Zi,Q) = E[E m€ « rKZ m ,e^(Z i ,e))y^'"F (0 ! ? m ) I Zt],p e Sf\ 



A°(i,p,Zt,Q) = D°(i,p,Zi, Q) I (Y^ p £ m6i ^o F° (0, t m ) ) 

(27) 

where A°(i,p,Zi,Q) represents the power to invest at date i for product p for a gas 
volume level consumed Q and a stochastic state vector z. As shown by [17|, the 
D°(i,p,z, Q) can be calculated during the dynamic programming recursion using : 

D (i,p,Zi,Q) = fi[l f+le5 , j ,4*(i + l.Zi+i,^i(Zi,f2))l^ i i , " 1 ^ (0,/+l) I Zi] 

+ £ E[D(i+l,p,Z i+1 ,Q + f(i,Z i ,Q))\Z i = z}. 
pes -'' 

The same procedure can be applied to derive a Delta estimation for the index com- 
modity component j > : 

' DJ(i,p,Zi,Q) = EE jDj{i,m,Z h Q)Fj(0,t m ) I Z>],P G J2/ 

P „• (28) 

AJ(i,p,Zt,Q) = DJ(i,p,Z„Q)/(Y? P Z mS£ ?jFJ(p,t m )) 



Using equation ( 15 i we get 

&(i,p,Z uQ )= £ £ F>(0,r M ) lr(t) - f(t)> ^- r(t) -* (t) E[^(^,/,Z,-,6)|Z,-] 

ffl£ #; * ie[r(t),r(*fi)[ mw 



Besides from equation ( 14 1 

E[D>(/,m,Z,Z„e) | Z /+1 ] = Eb*(/,Z /l Q;' i (Z /) !2))F^'" | Z, +1 ],form > i+l 

= E[q*(l,Z l ,Q*' i+1 (Z i+1 ,Q + q*(Z i ,Q)))Y^ m | Z i+1 ] 
= E[Dj{i,m,l,Z i+1 ,Q + q*(Zi,Q) | Z i+1 ] 

so that 

E[Dj(i,m,l,Zi,Q) | Z] - E[D>(/,m,Z,Z, +1 ,e + ? *(Z,e) | Zj] 

and 

D J (i,p,Zi,Q) = fi[l i+l6 ^(i,i+l ) Z,Zj,(2)) | Z,-] + 

£ fe[fl>(/+l,p,Z 1+ i,G + «*(Zi,fi))|Z>] (29) 

which can be easily calculated during backward recursion. 
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In first subsection we present the contract used for this study and we give the differ- 
ent numerical parameters used in calculation : 

1. the number of basis function used for each dimension for regression used to 
estimate conditional expectation J6], 

2. the number of particules used in an optimization part calculating the optimal 
value and the conditional hedge store in files, 

3. the number of particules used in a simulation part where the dynamic hedge is 
carried out. 

In a second subsection a study on the regression dimension is achieved and in a 
third subsection we test for the contract the sway of each component on the hedge's 
efficiency. At last we test the effect of the hedge frequency on the hedge efficiency. 



4.1 Contract description and general parameter 

We take the following example derived from a real Swing contract. We suppose 
that the Swing is contracted for year 2007 and that the day of valorisation is 1th of 
april 2006. The flexibility is such that each day normalized quantities are such that 

Qmax — 0.4, Q m in = 91,<2nwu = 146. 

Remark 4. With this kind of parameters the Q max constraint is ineffective. 

The arbitrage market taken for S° is the spot gas at Zeebruge hub. The index is 
composed of two components : 

• The first component is described by parameters (9,0,1) such that the average 
is taken on 9 months (m equivalent to 9 months), with no lag and the index 
component is changed every month. The first commodity is the brent. 

• The second component of the index is described by parameters (1,0,1) meaning 
that the average is taken during 1 month (m equivalent to 1 month) , no lag (Z = 0) 
and the index component 2 is changed every month. The second commodity is 
the month ahead gas on TTF market. 

The parameters a,, b\ have been adapted such that the Swing is at the money (the 
corresponding swap contract has a zero value) : ao = 2.2, a\ = 0.1757, ai = 0.2714, 
b x =45.16, b 2 =0. 

The model for each commodity are described by a two factor model (n — 2): 

• The Brent annual short term volatility is equal to 0.28, the annual short term 
mean reverting is equal to 0.1, and the annual long term volatility is equal to 0.1, 
the quotation is in dollar. 
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• As for the spot/future gas market, the short term annual volatility is equal to 1, 
the annual short term mean reverting equal to 80, while the long term volatility 
is equal to 0.1. The quotation is in Euro the domestic market. 

• Spread between the euro and dollar interest rate is 0, and the annual volatility 
between the two currency is 0.1 1. 

The daily forward curves are generated by Monte Carlo and at each day, depending 
on the market product available, this curves is averaged on each delivery period as- 
sociated to each product giving products values that will be used during hedging. 

The parameter used for the calculations are the following : 

• we use 80000 trajectories in optimization to calculate the Swing value and the 
hedge, 

• when regressing in dimension 1, 6 basis functions are taken. When regressing in 
dimension 2, 6 x 4 basis function are taken, while when regressing in dimension 
3, 6 x 4 x 1 basis functions are taken. 

• The global quantities Qi are discretized with a step of 0.4 and the local quantities 
with a step 0.1. 

Due to the calculation time and most of all due to the memory needed on computer, 

the Delta calculations have been parallelized by mpi (http://www.mcs.anl.gov/research/projects/mpi/) 

on a cluster during optimization and parallelization on scenarios has been achieved 

in simulation. A single optimization takes roughly 5 hours on 12 cores of a cluster 

and half an hour in simulation. 



4.2 Regression tests 

In the section we test three regressors : 

• for the first one, we approximate Z t by S^, 

• for the second one, we approximate Z, by (S^,I t ), 

• for the last one, we approximate Z, by /,,/,) where /, is the partial summation 
of index that will be used the following month 



Ii = Y^iijiS' 1 ', 




where 




Using the last approximation we take into account the fact that a new index is 
under construction for the next period. 
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In the last case, the I t and I t are very correlated and the regression procedure us- 
ing the Choleski method in the normal equation (see [6|) can fail when using many 
meshes in the last direction. That is the reason why we use only a linear approxima- 
tion in the last direction. 

Remark 5. The regression achieved is depending on the time step. If one regresses 
on Z t = (S® J t ) for a current step after t start , the regression is only achieved on S® 
before t start because /, does not exist. If one regresses on Z, = (S®,I t Jt) f° r a current 
step between t start and t enc i, Z, is restricted to (S^,I t ) before t start and (S®,I t ) after 

Results obtained in optimization and simulation are given in table [T] Standard de- 
viation of the portfolio with and without hedge is given. In this test, we hedge each 
component of the index and the exchange rate. As scheduled, the regression on 





Optimization 
value 


Simulation 
value 


Standard 
without hedge 


deviation 
with hedge 




86.2 


86.0 


219.3 


32.4 


z f = (s;u) 


97.7 


97.4 


204.1 


21.3 


Z, = {$,!„!,) 


99.6 


99.3 


202.9 


17.0 



Table 1 Tests on regressor 

Z, = Sf give bad results in term of value compared to the other regressors and in 
term of hedge efficiency. The regression on Z t = (S®,I t ) and Z t = (S^,I t Jt) have 
a far higher value indicating a better optimization than in the first case. The daily 
hedge is very efficient dividing the standard deviation by 10 for regressor 2 and even 
by 14 for regressor 3. 



4.3 Components to hedge 

In this subsection we choose to use a daily hedge calculated with regressor 3. On 
figure [3] we give for the 5 same simulations : 

• the hedging position in future October 2006 for Brent (so before the exercice have 
began) , 

• the hedging position in Brent future for january 2007 for the same simulations, 

• the hedging position for january 2007 in gas future, 

• the hedging position in foreign currency to hedge the product. 

In table [2] we give the standard deviation obtained in simulation taking different 
hedge : 

• for "Spot gas" we only hedge without taking into account the variation of the 
index components and exchange rate, 
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Delta on 5 scenarios 



Delta on a scera'los 
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(a) Hedging Brent position October 2006 
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(b) Hedging Brent position january 2007 
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(c) Hedging Gas position january 2007 
Fig. 3 Hedging position for different maturities 
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(d) Hedging exchange rate position 



• for "Index components without exchange rate" we only hedge against the varia- 
tion of the values of the commodities without hedging against gas spot variation 
and exchange rates variations, 

• for "Index components and exchange rate", we add to the previous hedge an 
hedge againt the exchange rate variations, 

• for "Exchange rate", we only hedge against the exchange rate variations, 

• for "Total hedge" we hedge against all the commodities and exchange rate in the 
contract. 

On this example, it shows that the hedge has to be done against the variations of 
the index components and that the hedge on gas spot variation is insufficient. It also 
shows that all the commodities and exchange rate has to be hedged to get an efficient 
dynamic hedging. On figure[4] we give de cash flow distribution obtained while not 
hedging, while hedging the index components, while hedging against all the product 
variations. 
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Hedge on 


Porfolio Standard deviation 


Spot gas 


196.61 


Index component without exchange rate 


77.20 


Index component and exchange rate 


74.03 


Exchange rate 


201.7 


Total hedge 


17.01 



Table 2 Tests on hedge 



No hedge 

With hedge 

With hedge on change and fuel in index 



0.030 



0.025 - 



0.020 



0.015 



0.010 



0.005 



0.000 




-400 



600 



Fig. 4 Test on the effectiveness of different hedges 



4.4 Frequency of the hedge 

It is well know that the hedging error in Black Scholes framework converges to zero 
at a rate proportional to the square root of hedging frequency lfT9lfT2ll . Because the 
index is evolving quite smoothly, we wonder if we could hedge the derivative less 
frequently as for the index components : it could help us to decrease the computa- 
tion time during optimization and would be interesting for practitioners because it 
would decrease the hedging cost due to illiquidity of the markets. In table[3] we give 
the standard deviation of the portfolio hedge supposing we hedge all the component 
but hedging the index with different frequencies. It shows that the effectiveness of 
the hedge decreases very quickly with the hedge frequency. In table |4j we give the 
standard deviation of the portfolio hedge supposing we hedge all the component but 
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Frequency of the index hedge 


Every day 


twice a week 


every week 


twice a month 


Standard deviation 


17.0 


42.9 


54.78 


57.5 



Table 3 Standard deviation of the hedged portfolio with different hedge frequencies for the index 
components 



hedging the index with different frequencies before the first exercice date and bal- 
ancing the hedging position of the index component every day after the first delivery 
day. Results in table[4]clearly shows that the hedge frequency can be lowered before 



Frequency of the index hedge be- 
fore t s t a rt 


every week 


twice a 
month 


every month 


every quarter 


Standard deviation 


18.6 


20.04 


20.79 


34.89 



Table 4 Standard deviation of the hedged portfolio with different hedge frequencies for the index 
components before the first exercice date and hedging the index components once a day after. 



the first delivery date without affecting too much the hedging efficiency. 



5 Conclusion 

We have derived an efficient methodology to hedge dynamically some Swing prod- 
uct on gas market. We have shown that the strategy calculated is very efficient for 
reducing the standard deviation of the portfolio for a realistic contract. Besides we 
have shown that on this contract the Swing value is mostly depending on the index 
component and that an efficient hedge has to deal with all the commodities involved 
in the index. At last, due to transaction cost it is always interesting to lower the hedge 
frequency and the last part of the study shows that is could be achieved before the 
first delivery date without losing too much the hedge efficiency. 
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